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Spontaneous formation of double bars in dark matter dominated 
galaxies 
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ABSTRACT 

Although nearly one-third of barred galaxies host an inner, secondary bar, the formation and 
evolution of double barred galaxies remain unclear We show here an example model of a 
galaxy, dominated by a live dark matter halo, in which double bars form naturally, without 
requiring gas, and we follow its evolution for a Hubble time. The inner bar in our model 
galaxy rotates almost as slowly as the outer bar, and it can reach up to half of its length. 
The route to the formation of a double bar may be different from that of a single strong bar 
Massive dark matter halo or dynamically hot stellar disc may play an important role in the 
formation of double bars and their subsequent evolution. 

Key words: galaxies: structure - galaxies: kinematics and dynamics - galaxies: spiral - 
galaxies: evolution - galaxies:halos 



1 INTRODUCTION 

A high fraction (> 60%) of disc galaxies in the local Universe are 
baiTed, including our Milky Way, and nearly 30% of baiTed galax- 
ies host an inner, secondary bar, nested inside the main bar, i.e., 
are do uble barred dErwin & SparkelBOOll : lErwinbOlU ILaine et al.l 
l2002h . The inner bars are likely to be old str uctures, as they are seen 
in the near-infrared dMulchaev et al.lll997l) and stellar population 
analysis giv es rather old age estimates ( Ide Lorenzo-Caceres et al.l 
I2OI2II2OI3I) . There is observational evidence that the two bars i n 
double barred systems rotate independently JCorsini et al.ll2003h . 
The observed common occurrence of double bars is not reflected in 
the present numerical models of bar formation. While a stellar disc 
with the Toomre parameter Q ~ 1 readily forms a single bar, which 
then grows within a Gyr to become a strong bar (such as the one in 
NGC 1300), long-term evolution of such a strong bar does not lead 
to the formation of a double bar on its own in purely stellar discs. 
Only a few N-body simulations have reported formation of double 
bars JRautiainen & Sal jl999l:ICurir et al.l2006tlDebattista & SherJ 
I2007D . mostly when special initial conditions ha ve bee n imposed. 

In 2D simulations of [Rautiainen & Sale ( Il999l) . with rigid 
bulge and halo, double bars form when the Toomre Q parameter is 
increased in the central parts of the disc, and they can sur vive for the 
Hubbl e time. In sets of cosmological A'^-body models bv lCurir et al.l 
(1200 6 )■ double bars form when the mass of the disc is lowered. 
lOeba ttista & Shen ( 2007) showed that the inner bar developing 
from a rapidly rotating bulge, or a pseudobulge, survives many 
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relative rotations of the bars. Such inner bar pulsates, and its pat- 
tern s£eed_£scU]^tesjnaMordwith_gredictions from orbital analy- 
sis JMacieiewski & Sparkdl 19971120001) . Another route to forming 
double bars relies on the pre sence of a dissipative (ga seous) com- 
ponent. In this scenario (e.g.. lFriedli & Martinetl ll993'). gas inflow 
in the large-scale bar stagnates in the inner kpc, leading to the for- 
mation of a disc there, which may become unstable and give rise to 
a smaller, secondary bar. However, in numerical realizations of this 
scenario, the inner bar last s no longer than a few relative rotat ions 
of the bars (see sect.5.2 of iMacieiewski & Athanassoulall2008l for 
a summary). 

In this letter, we report spontaneous formation of double bars 
in a dark matter dominated stellar disc without any gas. We have 
performed a suite of simulations of dark-matter-dominated galax- 
ies (Saha et al. in prep) in 3 dimensions that include a live halo. We 
noticed that in a few cases, structures resembling double bars form 
naturally in our simulations. In this paper, we present a model with 
such structure being most evident. It provides insight into factors 
decisive in formation of double bars, and explores self-consistent 
double bars of parameters markedly different from previous simu- 
lations. 



2 INITIAL GALAXY MODEL 

Equilibrium model of a galaxy is constru c ted us ing the self- 
consistent method of iKuiiken & Dubinskj 119951) . The initial 
galaxy model consists of a live disc, halo and a classical bulge. 
The disc has an exponentially declining surface density with a 




dSpringel et alJuOOlh with a tolerance parameter 9toi = 0.7, and 
the maximum value of the integration time step '^ 0.03, corre- 
sponding to 1.2 Myr. A total of 2.2 Million particles were used 
to represent the galaxy model with 1.05 Million each for the disk 
and halo and 0.1 Million for the bulge. 



Figure 1. Initial circular velocity. Red line denotes contribution from the 
disc, blue and green lines denote contributions from the bulge and the 
dark matter halo, respectively. Solid black line is the total circular veloc- 
ity. Length and velocity are given in internal units. See section|2]for scaling 
to physical units. 



scale-length Rd, scale-height hz, and mass Md- In internal units, 
where G=l, these parameters take the following values: Rd — 1, 
hz = 0.03 and Af^ — 1.58. The outer radius of the disc is trun- 
cated at G.ORd with a truncation width of O.SRd within which the 
stellar density smoothly drops to zero. T he live dark matter halo 
is modelled with a lowered Evans model JEvansll 19931) which has 
a constant density core. Such a cored halo is known to better rep- 
resent the observed high resol ution rotation curv es in low surface 
brightness (LSB) galaxies CKuzio de Narav et al..200q). The initial 
classical bulge is modelled with a King model (King 1 196a ). The 
mass of the dark halo is Mh = 20.43 and that of the classical 
bulge is A/f, — 0.153. For relevant details on model construction, 
the reader is referred to lSaha et al.l(l2010ll2012l) . The initial Toomre 
Q profile for the galaxy is such that Q rises to a high value beyond 
about 5Rd and the same happens at radius below ^ IRd- At2.5Rd, 
the Toomre Q reaches the minimum of Q = 2.55. 

In Fig. [T] we show the circular velocity curve for the galaxy 
model under consideration. The model galaxy is dark matter domi- 
nated right fro m the central region . This is a norm amongst most 
LSB galaxies jde Blok et alJuOOlr ). However, rotation curves in 
LSB galaxies usually show a slow rise, while the circular veloc- 
ity curve in our model rises sharply in the inner region. Such a 
sharp rise is seen in giant LSB galaxies which often contain a 
bulge-like component JBeiiersbergen et al.lll999l : lLelli et al.l2010r) . 
Our galaxy model has some resemblance to these giant LSBs, 
but unlike them it contains no gas; hence caution is advised if 
our results are to be used in studies of LSB galaxies. If we set 
the unit of length to Rd = 4.0 kpc and the circular velocity 
at R = 2.1Rd to 220 km s~^, then the units of time, mass 
and velocity are 42 Myr, 8.08 x 10^ Mq, and 93.2 km s"\ 
respectively. Dimensional values in the remainder of this paper 
are given in this standard scaling. Thus, in our standard scaling, 
the disc, bulge and halo masses are Md ~ 1.27 x 10^" Af©, 
Mb = 0.124 X 10^°Mq, and Mh = 1.65 x IO^^A/q, respec- 
tively. For any other mass unit Afo and length unit Rd, the time unit 
is 42 Myr(7?d/4kpc)^-^(A/o/8.08 x 10^ Mo)"^-^ and the veloc- 
ity unit is 93.2 km s"^(i?d/4kpc)~°'5(Afo/8.08 x 10^ Mq)" ^ 
Note that scale lengths of giant LSB discs are typically 10 kpc or 
more JBeiiersbergen et al.ll999l) , but for such scaling the simulated 
evolution time exceeds Hubble time. 

The simulation was performed using the Gadget code 



3 RESULTS 

The set of images presented in Fig.|2]represents density of the stel- 
lar component at different times throughout the run, projected onto 
the disc plane. The radial variation of the m = 2 Fourier com- 
ponent of the stellar density as a function of time is presented in 
Fig. [U In order to estimate the extent of the bars, and to measure 
their orientation, we fitted ellipses to the density field on a set of 
images such as in Fig. [2] and derived their ellipticity and the posi- 
tion angle (PA) in the same way as it is done for the observational 
data. In Fig. [4] we show the ellipticity and the PA as a function 
of the semi-major axis (SMA) obtained using the IRAF ELLIPSE 
fitting routine. If a bar is present, the PA of the major axis should 
be nearly constant over a range of sizes, with ellipticity reaching 
local maximum within this range. By matching the radial variation 
of the PA and the peak in the ellipticity, we assign an average value 
of the PA to each bar with an average error of '^ 10° . When the in- 
ner bar is not exactly perpendicular to the outer one, spiral features 
start from the end of the inner bar, making the measurement of the 
PA of the bar difficult. In this situation, we recheck our automated 
measurement of assigning a PA by eye. In me asuring the len gth 
of the bar, we follow the algorithm described by lErwinI fcOOSh for 
deriving Li,ar there. This is the upper limit for the length of a bar. 



3.1 Formation of two bars 

We follow the evolution of an initially axisymmetric stellar disc 
embedded in a dark matter halo, which gravitationally dominates 
the stellar component throughout the extent of the disc. Since the 
stellar disc is initially hot, it does not form a bar readily. There is 
no clear non-axisymmetric structure in the disc till t=48 in Fig. [2] 
which con^esponds to 2 Gyr, but a short open spiral can be no- 
ticed at t=72 (3 Gyr) in Fig.|2] The PA of the fitted ellipses at this 
time, shown in Fig. [4] increases almost monotonically with radius, 
hence the two bars are not well defined yet. However, there are 
two local maxima in ellipticity with values higher than 0.2, and 
the m=2 Fourier component in Fig. |3]also shows two maxima at 
t = 72. These two maxima, albeit with much lower amplitude, can 
be traced back to f = 48 at least, with the maximum correspond- 
ing to the outer bar forming first. Thus in our model two indepen- 
dent structural components are present from early stages of the run, 
which then develop into two bars. At later times, the two maxima 
in Fig.[3]correspond to the two bars. 

As the asymmetry grows in strength, the spiral transforms into 
two well defined bar-like structures that appear almost simultane- 
ously over time between t — 72 and 96 (3 and 4 Gyr). At i = 96, 
the two bars are nearly perpendicular to each other. Although there 
is still a spiral transition between the bars, the ellipse PA shown in 
Fig.|4]is roughly constant within the inner (up to 0.75 SMA length) 
and the outer bar (up to 1.5 SMA length). The outer bar has a clear 
boxy isophote. 



Formation of double bars 




Figure 2. Snapshots of face-on column density of all stars in our galaxy model. The color bai' represents density in logarithmic scale in internal units, which 
correspond to 505 Mq pc^^ in our standard scaling. The disc rotates counterclockwise. Snapshot time in internal units is shown in top-left comer of each 
panel. The first panel shows the initially axisymmetric disc that subsequently evolves into a double barred system. 




Figure 3. Time evolution of m = 2 Fourier component ( A2 ) of the stellar 
density distribution as a function of radius in polar coordinates aligned with 
the stellar disc. Color bar represents A2 normalized by the axisymmetric 
component Aq. The two bars are approximately confined within 2 scale 
lengths. 



3.2 Rotation of two bars 

After t — 48, the snapshots in Fig.|2]are shown every 24 time units 
(corresponding to 1 Gyr), which is approximately two rotation peri- 
ods of the outer bar (see below). The sequence clearly demonstrates 
that the inner bar rotates with respect to the outer bar. However, we 
find that the inner bar definitely is a very slowly rotating structure: 
it takes several rotations of the outer bar for the relative angle be- 
tween the bars to change considerably. During the period of 5 Gyr 
(between i = 96 and t = 216), the inner bar has rotated only once 
inside the outer bar, going from one state when the two bars are 
orthogonal to the next one. 

In order to quantify the rotation of the two bars, the PA of the 
outer bar, c&p, and of the inner bar, $5, were calculated every 0.3 
time units in the inertial frame. During every such interval, the PA 
of each bar increases by about 10°. In Fig. [S] we show how the 
phase difference of the two bars, $3 — $p, evolves with time. This 
difference increases monotonically, which means that the inner bar 
rotates faster than the outer bar. Past t — 190, the difference be- 
comes linear with time. In Fig.|6] we plot the pattern speed of each 
bar as a function of time. They are derived by fitting consecutive 
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Figure 4. Ellipticity and position angle of ellipses fitted at four instants 
(times shown in each panel) to the density distribution from Fig. 2 using 
IRAF ELLIPSE routine, as a function of the SMA length. At t = 72 the 
bars are not well defined yet, at t = 96 and t = 216 the bars are almost 
perpendicular to each other, and at t = 168 the bars are out of alignment. 



straight lines to <&p(f) and "I?s(i) data points over every period 
when each bar rotates by 360° in the inertial frame. These mea- 
surements are sufficiently accurate to imply that changes in pattern 
speed over time are real, though there is no clear regularity in these 
changes. The pattern speeds of the two bars are different only by 
~ 10% at the most. After t = 190, the rotation period of the outer 
bar is Tp = 2-k/Q,p ~ 12 time units. 

Although past t = 96 the two bars are separate entities, as 
demonstrated by nearly constant PA of fitted ellipses in Fig.|4l there 
is a spiral structure between them visible in Fig. |2] When the two 
bars are moving away from alignment at i = 136, towards be- 
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Figure 5. Phase difference (in degrees) between the two bars, 0g — ipp, as 
a function of time. 
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Figure 6. Pattern speed of the two bars as a function of time. Filled blue 
circles denote the primary bar and open black circles coiTespond to the sec- 
ondary bar. 



coming perpendicular at i = 216, this spiral is trailing. It turns 
to a leading spiral as soon as the bars get past the perpendicular 
arrangement at t = 216, and are on their way to become parallel 
again. Thus the spiral is always trailing when the bars are getting 
out of alignment, and leading when they are getting back to align- 
ment. The spiral structures are nearly absent when the two bars are 
perpendicular to each other. These characteristics are different from 
a spiral emerging at the ends of a bar that is driven by that bar, as 
in that case the spiral should always be trailing. The spiral in our 
model may be caused by the orient ation of orbits in the potentia l 
of the two bars, like in model02 of iMacieiewski & Smalll ( 120 lOh . 
when the loops (maps of orbits) form a trailing spiral when the bars 
are leaving the alignment, and a leading spiral when they are com- 
ing back to the alignment. When the bars are parallel or perpendic- 
ular, the loops are aligned, and therefore they do not form a spiral 
shape. The spiral structure may influence the dynamics of the two 
bars, and may indicate that the two bars are dynamically coupled, 
although in a different way than having resonances overlapping or 
pattern speeds commensurate. 



3.3 Evolution and dynamics of two bars 

The strength and size of the two bars increase with time, as can be 
seen in Fig. [3] and Fig. [4] Both of these quantities can be reliably 
measured when the two bars are orthogonal to each other, i.e. at 
t = 96 and 216. Between these two times, the length of the in- 




Figure 7. Resonances in the stellar disc. Black solid and dashed curves 
are axisymmetric approximations of f2 — k/2 and f2, respectively. Blue 
and green curves show Q — k/2 derived from tangential velocity along the 
outer and inner bar, respectively, when the bars are perpendicular (t = 96). 
The red and black hatched regions denote ranges of pattern speeds of the 
primary and the secondary bar, respectively, measured over the time period 
between t=90 and t=240. 



ner bar increases by two-fold: from 0.75Rd to 1.55Rd (see Fig.|4]l. 
In the same time period, the length of the outer bar grows mono- 
tonically from 1.5Rd to 2.1Rd. Thus the length ratio of the bars 
appears to increase from 0.5 at t=96 to 0.75 at i = 216. However, 
at relative bar positions other than orthogonal, the estimate of the 
length of the inner bar returns lower values, possibly because of 
the presence of a spiral structure connecting the bars, and then the 
length ratio remains close to 0.5. The amplitude of the peak in el- 
lipticity associated with the inner bar grows from 0.4 at t = 96 to 
almost 0.7 at t = 216 (see Fig.|4ll, i.e., by a factor of '^ 1.5. The 
ellipticity of the outer bar is lower than of the inner one: at i = 96 
and 216 it is about 0.3, although it reaches 0.5 at t = 168 and 
240. The increase of strength and size of the bars is moderated by 
the relative position of the bars: A2 within the outer bar is reduced 
when the bars are orthogonal, and A2 within the inner bar is re- 
duced when the bars are parallel (Fig.|3). The inner bar is growing 
stronger particularly after t — 144, as it grows in size and its A2 
increases. This is because of the combination of the secular and the 
periodic change caused by moving away from alignment. 

Having confirmed that the two bars are independent structures, 
one would like to know their dynamics: are they slow or fast bars, 
and what resonances they generate. Pattern speed of the outer bar 
decreases from about Sip — 0.58 at i = 136 to Sip — 0.50 at 
t — 237. Comparing these values with the azimuthal frequency 
curve in Fig.l?] we have the corresponding corotation radii at 3.6 
and 4.1Rd. On the other hand, the length of the outer bar increases 
in the same time interval from 1.66 to 2.25/?d. Thus the outer bar is 
slow, in the sense that it extends to only about 0.5 of its corotation 
radius, with this ratio increasing from 0.46 to 0.55. As the pattern 
speed of the inner bar is similar to that of the outer bar, the inner 
bar extends to even lower fraction of its corotation radius, with the 
ratio around 0.3. 

In order to determine the presence of the Inner Lindblad Res- 
onance (ILR), in Fig. [T] we plot the axisymmetric approximation 
to the fl — Hi/2 curve derived from the rotation velocity, accom- 
panied by the same curves derived from tangential velocity on the 
major axis of each bar, which relaxes the assumption of axial sym- 
metry. These curves do not differ much, which is expected in a 
model dominated by nearly spherical dark matter halo. Over the 
O — k/2 curves, we overplot the range of pattern speeds associated 
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with each bar throughout the run. Our measurements are consistent 
with either no ILR or a weak single ILR at around 1 — 1.2Rd, thus 
the inner bar cannot have its backbone built out of orbits related to 
the x2 orbits in the outer bar. Further work is needed to establish 
orbital support of double bars like the ones in the model presented 
here. The possible absence of an ILR makes the disc favourable to 
grow a bar through the swing amplification o f waves as it allows 
the feedback loop to complete l lToomrelll98lh . On the other hand, 
since our bars form slowly in an initially rather stable disc (high 
Toomre's Q), and since the y do not extend to th eir corotation radii, 
the mechanism proposed bv lLvnden-Belll ( 119791) may play a role in 
their formation. However, neither of these mechanisms anticipated 
formation of multiple bars. 

4 DISCUSSION AND CONCLUSIONS 



In this paper, we presented a model of a stellar disc, which spon- 
taneously forms two bars that is markedly different from systems 
simulated previously: the gravitational potential is dominated by 
the dark halo (Fig.[TJ, the inner bar is large (Fig.|2], and the angular 
velocities of the bars are almost equal (Fig.|6). 

All numerical simulations of double bars to date assume 
gravity dominated by stars in the region where the bars form. If 
Toomre's Q is low, then the outer bar forms rapidly, and an addi- 
tional process is needed to induce the formation of the inner bar. 
On the other hand, double bars can form spontaneously in pure 
N-body models when the disc is poorly coupled by its own self- 
gravity, which has to com pete with the gravity of the massive halo 
or with thermal motions. iRautiainen & Said 1 119991) were able to 
obtain double bars in their Model IV, in which they increased the 
Toomre parameter Q in the central parts of the disc to Q = 3, from 
Q — 1.5 in the otherwise identical Model I, which ret urned a sin- 

[ le ba r only. In sets of cosmological N-body models bv lCurir et al.l 
20061) ■ double bars form when the disc-to-halo mass ratio is small- 
est in each set. Only single bars form in more massive discs in 
those models. Our simulations presented in this paper confirm this 
trend, because they form double bars in the disc with high Q, which 
is dominated by dark matter halo. These findings indicate that the 
route to the fomiation of double bar may be different from that of 
a strong single bar, and the dark halo or hotter disc may play an 
important role. Bar formation in our simulation scaled to younger, 
smaller discs proceeds faster, but further work is needed to study 
evolution of such spontaneously formed double bars once the disc 
grows more massive. 

In the majority of numerical models, both purely stellar and 
with a gaseous component, the pattern sp eed of the irmer bar is sig- 
nifica ntly larger than that of the outer bar rFriedli & Martinelll993l: 
Rauti ainen & Salo 1999; Debattista & Shen 2007h iHeller et al.l 
2007D . Our model shows that an inner bar with the angular velocity 
similar to that of the outer bar is also possible. Since throughout 
the evolution of our model, the two bars can be in any relative ori- 
entation (F ig. [2](, the observed random orientation of the two bars 
JButa & C rocker 1993; Friedli & Martinet 1993) does not have to 
imply that pattern speeds of the bars differ significantly. 

The inner bar in our model is large - it is about half of the 
size of the outer bar for most part of the ru n. This is more than th e 
typical size ratio of the bars, being 0.12 dErwin & Sparkell2002h. 
but s ize ratios up to 0.4 have been observed (NGC 3358, lErwinI 
|2004) . An iimer bar supp orted by orbits inside the ILR of the outer 
bar cannot be too large (IMacieiewski & SparkdEOOOt) . but in our 
model the outer bar may have no ILR, hence other orbits, without 
such size constraint, must support the inner bar here. 



In summary, the model of double bars presented here indicates 
that (1) formation of double bars may proceed under different con- 
ditions and in a different way than the formation of a single strong 
bar - it may need dynamically hot stellar disc, possibly dominated 
by the dark halo; (2) inner bars as large as half of the length of the 
outer bar can last for a Hubble time or longer; (3) the difference 
between pattern speeds of the two bars can be minimal, yet the two 
bars can be observed in any relative orientation. 
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